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Abstract. We demonstrate the tunable enhancement of the zero phonon line of a 
single nitrogen-vacancy color center in diamond at cryogenic temperature. An open 
cavity fabricated using focused ion beam milling provides mode volumes as small as 1.24 
/im^. In-situ tuning of the cavity resonance is achieved with piezoelectric actuators. 
At optimal coupling to a TEMqo cavity mode the signal from individual zero phonon 
line transitions is enhanced by about a factor of 10 and the overall emission rate of 
the NV“ center is increased by 40% compared with that measured from the same 
center in the absence of cavity field confinement. This result is important for the 
realization of efficient spin-photon interfaces and scalable quantum computing using 
optically addressable solid state spin qubits. 


1. Introduction 

Coupling of fluorescence from nanoscale quantum systems to optical microcavities 
provides a means to control the emission process and can be an essential element of 
nanophotonic device applications. The negatively charged nitrogen-vacancy (NV“) 
defect in diamond is an example of a solid state system that has gained signihcant 
attention in recent years as a quantum spin register nca and nanoscale sensor mn 
due to its long spin coherence times and capacity for optical manipulation and readout. 
The recent demonstration of quantum error correction in an NV“ defect [B] provides 
a sound basis for using these systems in practical quantum information technologies. 
However a coherent spin-photon interface, necessary for many quantum technologies, 
can not be achieved if information is leaked via interaction with phonons, and so in 
the NV“ center only the ‘zero phonon line’ (ZPL) transition at 1.945 eV (637 nm) 
may be used. The low Debye Waller factor (DW ~ 0.04) of the center imposes 
a signihcant limitation, since most spontaneously emitted photons are useless. The 
enhancement of this transition and its efficient coupling to external optics are therefore 
important challenges to which microcavities are well suited. In particular is an essential 
requirement for generating large scale entangled states between spatially separated 


Tunable cavity coupling of the zero phonon line of a nitrogen-vacancy defect in diamond.2 

defects connected via photonic networks, for which proof of principle experiments have 
been achieved [8] but further development is hampered by low entanglement efficiency. 
Success in this endeavor will provide a route to the realization of scalable quantum 
computers based on optical networks of electronic and nuclear spins. [a HO] 

Resonant coupling of the NV“ ZPL to micro-ring resonators m and photonic 
crystal cavities [Ians] has been demonstrated to provide effective enhancements, but 
the monolithic structure of these cavities prevents positioning of the emitter at the 
heart of the cavity mode in-situ. Tuning of the cavity mode to the ZPL resonance 
can be achieved in these systems by progressive condensation of gas molecules onto 
the structure, thus increasing the mode volume and red-shifting the resonance, but it is 
difficult to optimise tuning using this procedure and it has limited scope for use in device 
applications. Open cavities [m [la [TB| [T7] provide a flexible approach to cavity-coupled 
devices that permit full in-situ alignment and tuning, combined with efficient coupling 
to external optics. Here we demonstrate control over the emission from a single NV“ 
center by coupling its ZPL to the resonant mode of an open microcavity. The open 
cavity geometry allows direct comparison between the emission properties of the same 
defect in and out of the cavity, thus providing unambiguous evidence of the effect of 
the cavity mode. We demonstrate clear tunable enhancement of the ZPL emission and 
reduction of the fluorescence lifetime of the defect in a controlled manner, and analyze 
our results in terms of the Purcell effect acting on the ZPL and the phonon sideband 
(PSB). Our work builds on recent demonstrations of room temperature coupling of NV“ 

(b) 1^ (d) 


I 


! 



Figure 1. Experimental configurations for characterisation of the same NV center 
in (a) the ‘out-of-cavity’ and (b) the ‘in-cavity’ measurements, (c) Scanning electron 
micrograph of a cross section through the Bragg coating of the concave mirror, (d) 
False color-scale plot of the electric field intensity distribution through a cross-section 
of the g = 4 TEMqo cavity mode calculated using FDTD modelling. The refractive 
index profile is included for clarity, with grey (purple) regions corresponding to the 
low (high) index layers of the Bragg mirrors. 
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defects to open cavities [IHlIiniEQl. 

2. Experimental Method 

The open microcavities are of a plano-concave design, the concave features produced 
by focused ion beam milling of a fused silica substrate [121 . These cavities support a 
Hermite-Gauss mode structure with TEMqo modes that can be effectively mode-matched 
to an external Gaussian beam. The radius of curvature of the concave mirror used here 
is 7.6 /im. The concave and planar dielectric mirrors have reflectivities of >99.99% and 
99.7% respectively at the design wavelength of 637 nm, and reflection bands extending 
from 550 nm to 720 nm. The planar mirror is terminated with a low index layer to 
provide a held anti-node at its surface. The shortest cavity length achieved here has a 
mirror spacing of 1.11 /xm providing an additional 3 held intensity maxima of the TEMqo 
mode between the mirror surfaces with A = 637 nm, as shown in hgure 1(d). We label 
this mode with the longitudinal index g = 4, or the set of indices (g, m, n) = (4, 0, 0). 

The NV“ defects are located in nanodiamonds of average diameter 100 nm produced 
using high pressure high temperature methods, which are spin cast onto the planar 
mirror. Registration marks on the planar mirror created using a focused ion beam allow 
individual defects to be identihed and characterised both in and out of the cavity. The 
two experimental geometries used in this report are displayed in hgure 1. The ‘out-of¬ 
cavity’ huorescence is measured in the absence of a concave mirror and with the planar 
mirror as a substrate behind the nanodiamonds (hg la). For the ‘in-cavity’ experiments 
the planar mirror is inverted and both excitation of the NV“ centers (at A = 532 nm) 
and collection of huorescence is carried out through the planar mirror. The microcavity 
assembly comprises a set of piezoelectric actuators that provide full control of the cavity 
length and relative position of the planar and concave mirrors at cryogenic temperature 
(hg lb). All measurements are carried out at 77K in a dry He exchange gas environment 
with the container supported in a liquid nitrogen bath cryostat IZH. Further details of 
the experimental apparatus can be found in section A of the supplementary information. 

3. Results and Discussion 

A low temperature huorescence image recorded in the out-of-cavity conhguration reveals 
well-dispersed nanodiamonds (hg 2a) with some instances of single centers recorded 
by the Hanbury Brown and Twiss method. The nanodiamond used in this study is 
circled, and its uncorrected photon intensity correlation data shown in hg 2b, revealing 
g^^^(O) = 0.38. Subtraction of background due to autohuorescence from the mirrors 
reduces this value to g*'^^(0) = 0.05 implying that ~ 97% of the huorescence is from a 
single emitter. Fluorescence spectroscopy of this NV“ center reveals a ZPL spectrum 
that is dominated by a linearly polarised doublet (htted peaks 2 and 3 in hgure 2d) 
indicating a high level of strain that lifts the degeneracy of the orthogonal and 
^Ey excited state dipoles of the defect (hg 2c). The line widths of the individual ZPL 
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components are approximately 0.4 nm, with Gaussian line shapes suggesting that they 
are dominated by inhomogeneous broadening due to spectral diffusion. Weak additional 
lines are seen, potentially due to an additional NV center in close proximity. The relative 
intensities of peaks 2 and 3 and the angle between their polarisations (fig 2e) indicate 
that the NV“ defect axis is oriented at 49° relative to the optical axis of the cavity, 
and the two orthogonal transition dipoles for lines 2 and 3 are at angles of 39° and 24° 
to the plane of the mirror respectively (see supplementary information section B). The 
ZPL measured has DW = 0.044, as is typical for NV“ centers. 


no" 




270 

Wavelength (nm) 


Figure 2. Characterisation of the NV“ center on the planar mirror, (a) Fluorescence 
image showing point-like emission from single NV“ centers, and registration lines 
generated by focused ion beam Ga’*' implantation that allow re-location of the same NV 
defect in the two experimental configurations, (b) Photon correlation data from the 
NV“ center measured, (c & d) Fluorescence spectra from the NV“ defect at T = 77K, 
and (e) polarisation plots of peak 2 (blue) and peak 3 (purple) of the zero phonon line 
doublet. 
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Figure 3. In-situ tuning of a cavity mode through the ZPL resonance, (a) Selected 
spectra (blue) recorded at five different cavity positions, with a simulation of each 
spectrum calculated using equation 5 (red), (b) 2D false color-scale plot of the 
measured spectrum including all of the tuning steps recorded. 

Figure 3 shows the measured fluorescence spectrum for the in-cavity configuration 
as a TEMqo cavity mode is tuned through resonance with the ZPL. Clear enhancement 
of the ZPL emission is observed, providing a fully saturated photon count rate of 15 
kc/s. This represents an experimental enhancement of the total ZPL signal by a factor 
of 2.5 compared with that recorded from the defect with the planar mirror alone (see 
supplementary information section C). 

Figure 4 shows a side-by-side comparison of the optical properties of the NV“ 
center out of the cavity, with that observed in the cavity at optimal tuning to the 
ZPL. Figure 4a shows the spectra over the full range of NV- emission, illustrating the 
extent to which the coupled ZPL dominates the measured fluorescence, a result of the 
cavity having no other modes in this range that couple efficiently into the objective 
lens. Figure 4b shows a comparison between the fluorescence decay characteristics. 
The lifetime of the out-of-cavity defect is measured as (30.8 ± 0.6) ns while that in the 
cavity is (22.1 ±0.4) ns, corresponding to a (39.5 ±0.7) % increase in the emission rate. 
The slight deviation from a single exponential decay in the in-cavity data at longer 
delay time is suspected to be due to spectral instability of the cavity mode leading to 
inhomogeneous broadening. Figure 4c shows the photon autocorrelation data, which 
reveal a reduced 5 '*'^^( 0 ) correlation value of 0.28 suggesting a slight improvement in the 
isolation of the single emitter when only the cavity-coupled ZPL is measured. 

We separate the analysis of our results into two parts - a semi-analytic treatment of 
the coupling of the ZPL to the TEMqo mode and a numerical treatment of the coupling 
of the PSB to other cavity modes present. We begin with the ZPL coupling, for which 
we express the wavelength-dependent enhancement of the emission of a given dipole by 
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Figure 4. Comparison of properties of the same NV“ center in the in-cavity and out- 
of-cavity geometries, (a) fluorescence spectrum (b) fluorescence decay, and (c) photon 
correlation histogram. In (a) and (c) the excitation power is set to Pgat as measured 
in intensity saturation measurements. 


a single cavity mode as: 
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factor between the emitting dipole /r and the cavity electric held E, Xcav is the cavity 
wavelength and Q is the cavity quality factor. As we are able to position the emitter at 
the electric held maximum we assume that E = E^ax so that = cos^ {9) where 9 is 
the angle between the dipole and the plane of the mirror. 

The fractional increase in the total emission rate arising from the coupling with the 
ZPL is given by 


F; 


ZPL 


= dXS^iX)F^iX) 


( 2 ) 


where are the branching factors of the two dipoles p, in the excited state and 5^ are 
their normalized free emission spectra collected over all directions. For peaks 2 and 3 of 
the ZPL doublet = 0.44 and 0.56 respectively (see supplementary information section 
B). Fzpl as dehned includes the Debye Waller factor through the relative weight of 
corresponding to the ZPL. We rewrite equation 2 as: 
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is the normalized spectrum emitted along the cavity axis, to which the spectrum shown 
in hgure 2(c) is a good approximation when appropriately scaled. The spectrum emitted 
into the cavity mode now reads 
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Figure 5. Comparison between (a) the measured in-cavity spectrum with the 
(5,0,0) mode optimally coupled to the ZPL, (b) the FDTD simulation of the Purcell 
enhancement, and (c) the semi-empirical calculation of the emitted power density 
spectrum. Modes are identified by their [q, m, n) for longitudinal index q and transverse 
indices to, n. 

which, taking SXcav = 0.7 nm, reproduces well the measured tuning spectra in £g 
3a. The same value of SXcav used in equation 3 gives Fzpl = 0.25. This quantity is 
equal to the ratio of the overall ZPL emission rate into the cavity mode with the total 
emission rate in free space. 

Based on the measured Debye Waller factor of 0.044 and the branching factor for 
peak 3 of 0.56, the cavity-induced enhancement of emission from peak 3 is found to be 
about a factor of 11. This compares with a value of 9.2 obtained using equation 1 with 
the effective Q factor Qe// = X/ {SXcavUy + SXemitter), in which SXcmitter = 0.4 nm, the 
fitted line width for peak 3. 
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Modification of the optical density of states experienced by the phonon sideband was 
determined using numerical Finite Difference Time Domain (FDTD) calculations that 
reflect the full dielectric environment of the NV“ center (see supplementary information 
section D for details). These calculations allow conhrmation of the cavity parameters 
by matching the simulated mode spectrum to that measured, and a direct prediction 
of the change in emission rate that will occur between the in-cavity and out-of-cavity 
experimental conhgurations. Figure 5 shows semi-logarithmic plots of three cavity- 
coupled spectra with the (5,0,0) mode tuned to the ZPL. This measured spectrum (£g. 
5a) is chosen as it reveals the positions of other modes coupling to the PSB and therefore 
provides a good reference by which to verify the geometrical parameters of the cavity. 
The relative Purcell factor (fig. 5b) and the resultant semi-empirical prediction for the 
emission power density spectrum of the NV“ center in the cavity (fig 5c) are also shown. 
The PSB emission couples to cavity modes with longitudinal index g = 4 and transverse 
indices 0, 2, and 4. The absence of observed coupling to modes with odd transverse 
indices suggests that the NV“ center is well positioned on the cavity axis of symmetry 
where the electric held intensities of these modes drops to zero. 

Integrating the emitted power density spectrum between 640 nm and 740 nm reveals 
the Purcell enhancement of the NV“ center PSB emission to be = 0.93 (ie, a 

7% suppression of emission relative to the out-of-cavity geometry). The total change in 
emission rate for the NV“ center is thus predicted to be Fzpl+Fpsb = 1.14 significantly 
lower than the measured value of 1.395. 

Figure 6a shows the calculated Purcell enhancement factors of the ZPL and PSB 
for different TEMqo modes tuned into optimal resonance with the ZPL. The ZPL data 
points correspond to the enhancement of the emission rate corresponding to the entire 
ZPL: an estimate of the enhancement of peak 3 alone can be obtained by multiplying 
these values by l/n^ = 1.8. Figure 6b reveals that the theory above underestimates 
the experimentally observed lifetime changes for each of these cavity lengths. The 
total emission rate is predicted to remain approximately the same for g = 4, 5 and 6 
because the increase in ZPL emission as the cavity length is reduced is compensated by 
a reduction in PSB emission. 

We attribute the difference between the predicted and measured behavior to the 
presence of inhomogeneous broadening, both in the NV“ center ZPL due to spectral 
drift and in the cavity linewidth due to mechanical instability in our apparatus. Under 
such conditions the fastest component of the measured emission decay curve corresponds 
to near-resonant alignment and will primarily reflect the homogeneous line widths. We 
introduce inhomogeneous broadening of the cavity mode into our analytic calculations 
above by convoluting the mode line width with a Gaussian function g{Xcav)- The change 
in the NV“ emission rate due to the ZPL coupling is then given by (supplementary 
information section E); 


J dXcav g(^Xcav^(^7 FzPL{Xcav'))FzPL{Xcav) 

^ fdX^a. 9{Xca.)FzPL{Xca.) ^ ^ 

A cavity line width of SX = 0.2 nm, corresponding to the value measured 
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Figure 6. Calculated and measured parameters as functions of the longitudinal mode 
index q. (a) Calculated Purcell enhancement factors for the zero phonon line (open 
stars) and the phonon sideband (closed triangles), (b) Comparison of the measured 
emission rate enhancement (solid squares) with the calculated values assuming a cavity 
line width of 0.7 nm and no inhomogeneous broadening (open circles) and with a cavity 
line width of 0.2 nm and inhomogeneous broadening of 0.5 nm (open triangles). 


by transmission spectroscopy in a nominally identical cavity, combined with an 
inhomogeneons broadening of 0.5 nm, gives FzpL,inhom = 0.364, whilst leaving the 
detnning spectra in hgnre 3 relatively nnchanged. The overall changes in the NV 
emission rates are plotted in hgnre 6b, and are seen to agree more closely with, althongh 
still consistently nnderestimate, the measnred valnes. We attribnte the remaining 
discrepancy between the measnred and modelled lifetimes to inhomogeneons broadening 
of the ZPL, which is more difhcnlt to qnantify and which the semi-empirical calcnlation 
method described above can not easily accommodate. A simple indication of the 
potential effect to the lifetimes in this experiment can be obtained from eqnation 1 
however, in which a ZPL homogeneons line width of order 0.1 nm, consistent with 
valnes measnred at this temperatnre in bnlk materials [23], and a cavity line width 
of 0.2 nm, give = 33.6 when resonantly tnned to peak 3. The resnltant emission 
rate increase is Ftot = 1.71, snggesting that the measnred rate increase can indeed be 
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accounted for by a combination of cavity and ZPL inhomogeneities. For completeness 
we note that an unknown parameter for our NV“ center is its quantum efficiency (QE), 
which the above treatment assumes to be unity. Recent reports have demonstrated 
that nanodiamonds of the size used here can have QE as small as 0.3 [25]. The effect 
of reduced QE would be to require higher still Purcell factors to achieve the measured 
lifetime change, since the enhancement acts only on the radiative term. 

In conclusion we have shown the controlled coupling of the ZPL of a single NV“ 
center in nanodiamond to an open cavity at cryogenic temperatures. This cavity system 
shows signihcant potential for interfacing the NV“ center with photonic networks and 
performing quantum operations and measurements. The degree of enhancement is 
currently limited by the line width of the NV“ centers in the nanodiamond, and by 
the cavity line width as determined by scattering and instability of the low temperature 
cavity assembly. Improvements in these areas are readily achievable and are expected 
to produce much larger enhancements than those reported here. Open cavity quality 
factors exceeding 10® have been demonstrated irzi, and modest further reductions in 
cavity mode volume are clearly also possible. Single NV“ defects implanted into bulk 
diamond at depths of 100 nm can offer ZPL line widths as narrow as 27 MHz, |2B] 
and electron spin coherence times > 100/is. Resonant ZPL coupling of NV“ centers 
situated in diamond membranes to open cavities can thus result in enhancement of 
the emission rate into the ZPL by a factor of > 10^, leading to an effective Debye 
Waller factors approaching unity and indistinguishable photons with lifetimes of a few 
hundred picoseconds. Such projections suggest that the experimental conhguration 
demonstrated here is an attractive route towards an efficient spin/photon interface and 
to the construction of scalable quantum processors. 

This work was funded by the European Union Seventh Framework Programme 
(FP7/2007-2013) under grant agreement no 61807. SJ acknowledges support from the 
United Kingdom Engineering and Physical Sciences Research Council. 
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Supplementary Information 

A) Experimental Methods 

Cavities are milled into fused silica substrates (UQG Optics) using a focused ion beam 
(FEI FIB200). High-reflectivity mirror stacks are then deposited on the substrates. The 
planar mirror is coated in-house using alternating A/4n layers of Si02/Ti02 to achieve 
99.7% reflectivity at A = 637nm. The planar mirror is terminated with the low index 
materials such that there is an anti-node of the electric held at the surface for emitter- 
cavity coupling. The concave mirror was coated at LaserOptik GmbH, with 20 pairs of 
Si02/Ta205 to provide R/,99.99%, such that optical coupling is preferentially through 
the planar mirror. Nanodiamond solutions (lOOnm - O.lmg/ml) are spin coated onto the 
planar mirror substrates prior to cavity coupling. Optical characterisation is carried out 
with a home-built beam-scanning confocal microscope employing a fast-steering mirror 
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(Newport FSM300), outlined in figure 532nm CW excitation is used for imaging and 
spectroscopy, with 532nm pulsed excitation available for PL lifetime experiments (Teem 
Photonics-SNG-20F-lSO, repetition rate = 20KHz). Excitation is coupled through a 
single mode fiber (Thorlabs SM460HP) with polarisation control (Thorlabs FPC030). 
The cavity apparatus is situated in a dry He-exchange gas environment (50 mbar), 
immersed in a liquid nitrogen bath cryostat with optical access. The apparatus consists 
of a custom sample stage, allowing independent piezo actuation of the mirror substrates. 
The planar mirror substrate has all the translational degrees of freedom (Attocube: 2 x 
ANPxlOO, 1 X ANPzlOO), whilst the cavity substrate can only move vertically (Attocube 
ANPz30). All stepper motors are driven with an Attocube ANC300 control module. A 
low temperature compatible achromatic objective (Attocube ASWDO x50 0.82NA) is 
used for optical excitation and collection. Spectral filters can be placed in the collection 
path. Cavity coupled HBT, Lifetime and power saturation measurements are taken in 
the 633-647nm spectral window using an additional band pass filter (Semrock FFOl- 
640/14-25). Fluorescence is coupled into single photon avalanche detectors (Perkin- 
Elmer SPCM-AQRH) via a single mode hber (Thorlabs SM600). A 1x2 50:50 hber 
splitter (Thorlabs FCMM50-50A) with second single photon detector are used for 
HBT measurements. An Edinburgh Instruments TCC900 card provides the correlation 
electronics for the HBT and lifetime measurements. The output fiber can also be coupled 
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Figure 7. Optical set-up for low temperature cavity experiments 
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to a spectrograph (Acton 2500i), equipped with a liquid nitrogen cooled CCD (Princeton 
Specl0:100B). 


B) Calculation of the ZPL dipole orientations 

The orientations of the dipoles for the two strain-split ZPL transitions are calculated 
by determining the polar angle of the NV axis relative to the optical axis {6) and the 
rotation of the dipoles about the NV axis (/9) which will lead to the projections on the 
observation plane. 

A thermal distribution between the populations of the two excited states comprising 
the ZPL doublet is also assumed, consistent with the findings of Fu et al [23]. The 
energy splitting of peaks 2 and 3 is 1.5 meV, so that the population ratio for the levels 
responsible for peaks 2 and 3 at 77K is 0.8:1. The measured intensity ratio between 
peaks 2 and 3 due to the projection of this thermally distributed dipole pair onto the 
measurement plane is 0.58:1, so the equivalent projection of a circle formed by two 
perpendicular dipoles of equal strength would result in a ratio oi R = 0.73. 

The polar angle 9 is obtained from the sum of the polar intensity distributions of 
peaks 2 and 3, whereby the extrema have the following dependence on 9. 

hmin hmaxCOS (0) (7) 


The combined rotation matrix for the axial, then polar, rotation is 


A 


^ cos(/3) 

cos{9)sin{(3) 
^ sin{9)sin{(3) 


—sin{l3) 

cos{9)cos{(3) 

sin{9)cos{(3) 


0 ^ 
—sin{9) 
cos{9) j 


( 8 ) 


Applying this rotation to unit vectors X & Y, and projecting the resultant vectors 
onto the measurement plane, gives 

= 1 — sin^{9)sin^{(3) Y"^ = 1 — sin^{9)cos^{(3) 


Where the square of the dipole vector have been taken to obtain the intensity. 
Knowing 9 and the ratio R between these intensities, allows jd to be determined. Finally, 
the angles cfx' & (fv between the rotated dipole vectors and the observation plane are 
found by taking the dot product of the rotated dipole vectors and their projections, 
leading to 

(fx' = arccos f — sin'^{9)sin^{j3)\ cfy' = arccos (^Jl — sin^{9)cos'^{(3) 


C) Comparison of excitation conditions for the two experimental geometries 

For comparison of the emission intensities in the out-of-cavity and in-cavity geometries 
it is necessary to establish equivalent excitation conditions in the two cases. To do this 
we measured the dependence of the emission intensity I on excitation power P to record 
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the emission saturation curve of the color center in each case. These data are shown in 
hgure 8. We then htted these curves to the saturation function 

j _ I sat P 

Psat + P 

The htting parameters Igat and Psat are the fully saturated photon count rate 
and the characteristic saturation power for excitation respectively. In the out-of-cavity 
geometry we hnd that Igat = 154 kc/s for the whole NV centre (6.78 kc/s for the ZPL 
only) and Pgat = 1-02 mW, whilst that in the q = A cavity (data in hgures 3 and 4) 
yields Igat = 15.1 kc/s for the ZPL only and Psat = 1-89 mW. The spectra shown in 
hgure 4 use an excitation power of Psat in each case. 


a) Out of Cavity 


B 15000 
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Figure 8. Power Saturation measurement into the cavity coupled ZPL with mode 
number q=4. 


D) Finite Difference Time Domain simulations 

We used Lumerical Solutions FDTD software to perform numerical simulations of a 
dipole emitting into the cavity. Firstly, by reproducing the measured mode structure 
using these simulations we conhrm the cavity geometry. The concave mirror is 
modelled as vertically stacked hlms of A/4 vertical height, based on cross-sectional SEM 
measurements such as the one shown in hgure 1(a). Refractive indices of the Si02 and 
Ta 205 layers are 1.52 and 2.10 respectively. The concave mirror is situated as having 15 
Bragg pairs rather than the 20 in the real device, as this allows faster calculations with 
negligible ahect on the results. The dipolar source is positioned on the axis of symmetry, 
20 nm from the planar mirror surface, and within a dielectric sphere of diameter 100 nm 
and refractive index 2.4 in contact with the planar mirror to simulate the nanodiamond. 
Separate calculations are performed with dipoles at angles of 18° and 39° to the planar 
mirror representing peaks 3 and 2 respectively. We conhrmed that the results of the 
calculation are insensitive to the position of the source within the nanodiamond and on 
the exact size of the nanodiamond. The hnite element Yee cells have a minimum size 
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of 5 nm to accommodate the nanodiamond and the contours of the curved mirror. The 
electromagnetic held was allowed to propagate and decay for 10 ps in 5 fs increments, 
corresponding to a simulation time of about 12 hours. 

Mode volumes were calculated both using the simple analytic formula for a Gaussian 
beam and numerical integration of the simulated FDTD held intensity. The analytic 

expression used is Knode = where zr = — 1^ is the Rayleigh range, L is the 

optical length of the cavity and (5 is the radius of curvature of the concave mirror. The 
integration used is 

JRe[e]\EfdV 
(Re[e] \Ef) 

\ / dipole 


E) Calculation of inhomogeneous broadening 


For an inhomogeneously broadened cavity mode, the total intensity is given by an 
average over the diherent cavity positions 

I{t) = J dX,a. 9{Xcav)A{X,av)e-^‘'^-'>^ (9) 

where g{Xcav) is the inhomogeneous distribution of the cavity position, taken as a 
normalized Gaussian. For a given cavity position, the fraction of light emitted in the 
cavity mode is proportional to the measured intensity integrated over a full decay event 

;3(Ao„.) = , a /1(A„.)/7(A„„) (10) 

i + rzPLy^cav) 

The decay rate is then found by diherentiating equation 9: 


so that 


__dl 

'tinhorn — 

Idt t=0 


inhom 


'tinhorn 

7o 


_ fdX cav 9{^ cav )7^(A cav )P{X 

cav) 

JdX 

cav cav )7(A cav )P{X 

cav) 

f dXcav ^7(Acati)(l T ^ZPriXcav)^ fdi^Xoayj 
J dXoav fi^(Acati)(l T (Acan))/3(Acan) 


( 11 ) 


( 12 ) 


whereby substituting back in for /3{Xcav) from equation 10 gives 
^ J dXcav fi^(Acati)(l T RzPL(A cati))R zPL(A can) 


inhom 


JdX 

cav ^(A cav )Ezpl{X cav) 


(13) 



